Realistic description of electron-energy loss spectroscopy for one-Dimensional Sr 2 Cu03 
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We investigate the electron-energy loss spectrum of one-dimensional undoped CUO3 chains within 
an extended multi-band Hubbard model and an extended one-band Hubbard model, using the 
standard Lanczos algorithm. Short-range intersite Coulomb interactions are explicitly included in 
these models, and long-range interactions are treated in random-phase approximation. The results 
for the multi-band model with standard parameter values agree very well with experimental spectra 
of Sr2Cu03. In particular, the width of the main structure is correctly reproduced for all values of 
momentum transfer. We find no evidence for enhanced intersite interactions in Sr2Cu03. 

PACS numbers: 71.27.+a, 71.45. Gm, 71.10.Fd 



One-dimensional systems are easy to conceive in theory 
but hard to find in nature, and their experimental realiza- 
tion is restricted to few materials. These, include meso- 
scopic systems like single- wall nanotubead'EI or chains of 
metal atoms and macroscopic systems with a strong 
anisotropy in one spatial direction. Among the latter, 
Sr2Cu03 has been in the focus of recent research. It 
contains separated chains of corner-sharing CUO4 pla- 
quettes, and is related to high-temperature supercon- 
ducting compounds of higher dimensionality. Gener- 
ally, the electronic properties of Sr2Cu03 are dominated 
by strong correlations of the valence holes. The low- 
dimensional character of magnetic excitations in this ma- 
terial manifests itself in magnetic susceptibility measure- 
ments that have been successfully described in termSppLa 
one-dimensional spin-i Heisenberg antiferromagnet.tJu 

Charge excitations in S^CuOa have recently been 
investigated!!] by means of electron energy-loss spec- 
troscopy (EELS). The experimental spectra are shown 
in the right panel of Fig. [|. In the following we will only 
discuss the spectral region below 4eV energy loss. Ex- 
citations at higher energies probably involve Sr orbitals, 
which are not included in models for the Cu-0 structured 
In the low energy region the experimental data show a 
broad dominant low-energy structure at 2.4 eV for mo- 
mentum transfer q = O.lA -1 , which shifts to 3.2eV for 
q = 0.8A -1 [see Fig. The behavior of this struc- 
ture as a function of momentum transfer is rather un- 
usual: with increasing momentum transfer up to 0.4A -1 
the width of the structure decreases but increases again 
for q > 0.4A -1 . These spectra have been interpreted us- 
ing an extended one-band JJubbard modeljj an effective 
two-baiid Hubbard model,£3 and a multi-band Hubbard 
model.Eil Although the main difference between these 
models is just the elimination of the oxygen degrees of 
freedom, the results have been discussed controversially. 
In the one-band model the behavior of the low-energy 
feature has been interpreted as the transfer of spectral 
weight from a continuum of excitations to an exciton, 
formed due io a strongly enhanced intersite Coulomb re- 
pulsion vJj'tHl Excitonic features have also been discussed 
in the strong coupling limit of an effective two-band Hub- 



bard model. EHl However, the used coupling strengths are 
not experimentally relevant. Therefore, no direct com- 
parison to experimental data was possible in Ref. Eo[ In 
contrast, in the multi-band model the dispersion of the 
low-energy feature has been explained in terms of a shift 
from a rather delocalized Zhang-Rice singlet-like exci- 
tation to a more localized oneO However, no intersite 
Coulomb repulsion was included in the Hamiltonian. 
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FIG. 1. Comparison of experimental data for S^CuOs 
(right panel), taken from Ref. |^, and the results of the exact 
diagonalization (left panel). The parameter set is Ud = 8.8eV, 
A = 3.0eV, V pd = 1.2eV, V dd = OeV, t pd = 1.3eV, and 
t pp = 0.65eV. The theoretical line spectra have been con- 
voluted with a Gaussian function of width 0.35eV. 

Up to now, all theoretical approaches failed to cor- 
rectly describe the discussed decrease and increase in 
spectral width of the low-energy feature as a function 
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of increasing momentum transfer: For small momentum 
transfer, the one-band model overestimates the experi- 
mentally observed width by a factor of about twola In 
addition, the broadening for large momentum transfer 
is to small. The analytical approach to the multi-band 
model, on the other hand, underestimates the broaden- 
ing due to the neglect of far reaching excitations that are 
important at small momentum transfer O 

In this paper, we show that the multi-band model pro- 
vides a realistic description of the EELS spectrum for 
Sr2Cu03, and we observe the correct spectral form for all 
values of momentum transfer. Furthermore, it is found 
that the main effect of the intersite Coulomb repulsion is 
to lead to an energy shift of the EELS spectra. Finally, 
we discuss the relations of our results to the loss function 
of the one-band model. 

We investigate the dielectric response of a one- 
dimensional extended multi-band Hubbard Hamiltonian 
at half-filling, i.e. a chain of corner-sharing CUO4 pla- 
quettes with one hole pex .Cu site. In the hole picture 
this Hamiltonian readsE£lLd 

ja i 

+V pd ]T rfnf + V dd ]T nf4 

<ij> <W> 



<ij>a 

+t pp ]T 4>^p] a P3' 
<jj'><? 



(1) 



The operators d icr iPj a ) create a hole with spin a in the 
i-th Cu 3d orbital (j-th O 2p orbital), and nf a (n?) are 
the corresponding number operators, with nf — nf^ +nf^ . 
The first four terms in Eq. ([!]) are the atomic part of the 
Hamiltonian, with the charge-transfer energy A, the Cu 
on-site Coulomb repulsion U d , the Cu-0 intersite repul- 
sion V pd , and the Cu-Cu intersite repulsion V dd . The last 
two terms in Eq. ([!]) describe the hybridization of Cu 3d 
and O 2p orbitals (hopping strength t pd ) and of O 2p 
orbitals (hopping amplitude t pp ). 4> % pd an d ^pp §i vc the 
correct sign for the hopping processes, and (ij) denotes 
the summation over nearest neighbor pairs. Hamiltonian 
(|l|) takes account of both in-chain and out of chain oxy- 
gen sites. Notice that no perturbative approximations 
are made, so that parameter values can be chosen in the 
experimentally relevant range. 

The dynamical density-density correlation function 
is directly proportional to the loss function in EELS 
experiments.ll3 By including the long-range Coulomb 
interaction in—, the model within a random-phase 
approximation!!!} (RPA) one finds for the loss function 



L(u>, q) = Im 



-1 



! + «qX»(w,q) 



(2) 



X°(w,q) = jrj dt 



e-'<(% q (i),p- q ]|0> 



(3) 



is the response function at zero temperature of the 
short-range interaction model ([|). X° p depends on the 
energy loss u) and momentum transfer q. |0) is the 
ground state, p q denotes the Fourier transform of m, and 
Vq = e 2 N/(eoe r vq 2 ) is the long-range Coulomb interac- 
tion with unit cell volume v. N is the number of electrons 
per unit cell, and e is the permittivity. The real part e r 
of the dielectric function can be obtained from the expeiu 
iment. In the case of S^CuOa it was found to be e r = 8. El 
In the following we_evaluate Eq. (^) using the standard 
Lanczos algorithm!!!! which is limited to small clusters. 
The theoretical line spectra are convoluted with a Gaus- 
sian function of width 0.35eV, to allow a comparison with 
the experiment. 
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FIG. 2. Finite-size effects in the loss function of model (|^) 
for clusters with six plaquettes (full lines) and five plaquettes 
(dashed line), with open boundary conditions. Parameters as 
for Fig. 0. 

First we check if our results are sufficiently converged 
with respect to system size. In Fig. || we compare the 
loss function of clusters with five (dashed lines) and six 
plaquettes (full lines). In both cases open boundary con- 
ditions are chosen. One has to make sure that holes on 
the edges of the cluster are still embedded in the local 
Coulomb potential that results from a state with occu- 
pied Cu sites. For this purpose, O (Cu) sites on the edge 
of the cluster are assigned an additional on-site energy 
due to V pd (V dd ). As can be seen from Fig. there are 
only small finite-size effects. Thus we conclude that the 
cluster with six plaquettes is large enough to obtain reli- 
able results. 
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In Fig. |l]the calculated loss function is compared to the 
experimental spectra from Ref. The parameters in the 
model Hamiltonian are chosen as follows: U d — 8.8eV, 
V pd = 1.2eV, V dd = OeV, t pd = 1.3eVjnd t pp = 0.65eV 
are kept constant at typical valuesJl3~E3 The value of 
A = 3. OeV has been adjusted to obtain correct peak po- 
sitions. This means that we use only one free parameter. 
As compared to the standard value 3.5eV for A,Ej Ej the 
smaller A is in agreement with theoretical -ajialysis of 
x-ray photoemission spectra for Sr 2 Cu03.at3E3 Notice 
that a small value of A means that the system is not 
in the strong coupling limit as was assumed in Ref. [l^. 
The calculated loss function consists of a dominant struc- 
ture at 2.5eV for q = O.lA -1 , which shifts to 3.2eV for 
q = 0.8A -1 . Besides, a second excitation is observed 
at 5.5eV. In agreement with the experimental observa- 
tion, with increasing momentum transfer the main struc- 
ture shifts to higher energies and first decreases in width. 
For q > 0.4A -1 the spectral width increases again. The 
main structure results from excitations in which a hole 
leaves its paginal plaquette to form Zhang- Rice singlet- 
like stateso with neighboring holes. With increasing mo- 
mentum transfer q the spectral weight shifts from ex- 
tended to more localized excitations.Eil 
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FIG. 3. Influence of V p d, A, and Vdd. on the loss func- 
tions with momentum transfer q = 0.3A -1 (left side) and 
q = 0.7A -1 (right side); other parameters as for Fig.nl. In 
(a) A = 3. OeV and Vdd are kept constant and V p d is varied 
(full line: V p d = 0, dashed line: V p d ~ l.OeV), in (b) A is 
varied (full line: A = 3. OeV, dashed line: A = 4.0eV) for 
V p d = Vdd = 0, and in (c) A = 3. OeV and V p d = are con- 
stant (full line: Vdd = 0, dashed line: Vdd = 0.5eV). The 
loss functions with q = 0.3A -1 and q = 0.7A" 1 are scaled 
independently of each other. 

Next we discuss the dependence of the calculated spec- 
tra on the different parameters in Hamiltonian (|l|). In 
agreement with an analysisc3 of the optical conductiv- 
ity for Sr 2 Cu03 it is found that the main influence of A 
and the intersite Coulomb repulsion V pd is a shift of the 
excitation energy of the main structure. This is shown 



in Fig. H for V pd [see panels (ai) and (ag)] and A [see 
panels (bi) and (b 2 )]. Note that we observe nearly the 
same results for A = 3. OeV, Vpd = l.OeV [dashed lines in 
panels (ai) and (a 2 ) of Fig. |j and A = 4. OeV, V pd = 
[dashed lines in panels (bi) and (b 2 ) of Fig. Hence, 
only the sum of both parameters A + V pd is relevant 
for the spectra. Therefore, it is possible to obtain good 
agreement between experimental and theoretical spectra 
with or without intersite Coulomb interaction V pd . This 
implies that the mechanism of excitations is not driven 
by a strong intersite interaction V pd . Furthermore, this 
observation also explains why a fit of a multi-band model 
with V pd — to the experimental data has led to a larger 
value of A in Ref. 
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FIG. 4. Influence of the Coulomb interactions U and V 
of the extended one-band Hubbard model on the loss func- 
tion with momentum transfer q = 0.3A -1 (left side) and 
q = 0.7A -1 (right side); parameter values are chosen accord- 
ing to Ref. pi The hopping strength is t = 0.55eV. In (a) 
V — 1.3eV is kept constant (full line: U = 4.2eV, dashed line: 
U = 4.7eV), and in (b) V is varied (full line: V = 1.3eV, 
dashed line: V = 0.8eV) for U = 4.2eV. The theoretical 
line spectra have been convoluted with a Gaussian function 
of width 0.35eV. The loss functions with q = 0.3A -1 and 
q = 0.7A" 1 are scaled independently of each other. 

In contrast to A and V pd , the Cu-Cu intersite repulsion 
V dd does not shift the complete EELS spectra but rather 
transfers spectral weight to excitations with smaller en- 
ergy loss [see panels (ci) and (c 2 ) of Fig. ||. A compari- 
son of panels (ci) and (c 2 ) shows that this effect is larger 
near the zone boundary at q = O.SA^ 1 . Probably, this 
behavior can be connected with a formation of an exciton 
state as discussed in Ref. |9|,|l0|- On ther«ther hand, the 
choice V dd = is a good approximationEa for the multi- 
band model (Q) since the distance between neighboring 
Cu sites is relatively large. Therefore, a possible exciton 
formation seems not to be relevant for the interpretation 
of the experimental spectra if one uses the multi-band 
model. 

Finally, we want to discuss the relation of our results 
for the multi-band model to the loss function (^) of the 
one-dimensional extended Hubbard model 
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and SFB 463. The calculations were performed on the 
Origin 2000 at Technische Universitat Dresden. 



(4) 

which only considers effective Cu 3d orbitals. In Eq. (|4|), 
t is the hopping strength, U denotes the on-site Coulomb 
repulsion, and V is the intersite interaction. Note that 
the mapping of the multi-band model (|l|) onto the one- 
band model (Q) is problematic for the model parameters 
used above because the system is not in the strong cou- 
pling limit. Therefore, we may compare the loss function 
of both models only qualitatively. In the following, we 
compute the loss function (Q) of the one-band model (ji|) 
using a cluster with twelve sites and periodic boundary 
conditions. If one reduces the multi-band model ([!]) to 
a one-band model the charge-transfer gap A is-seplaced 
by the Hubbard gap U of the effective model.E3 There- 
fore, in analogy to the influence of A in the multi-band 
model, increasing U shifts the spectra to higher energy 
loss [see panels (ai) and (a 2 )of Fig. In panels (bi) 
and (b 2 ) of Fig. || the loss functions with q = 0.3A -1 
and q = 0.7A -1 are shown for V = 1.3eV (full lines) 
and V = 0.8eV (dashed lines) where U = 4.2eV and 
t = 0.55eV are kept constant at the values from Ref. ||. 
In analogy to the influence of the intersite Coulomb re- 
pulsion Vdd in the multi-band model discussed above, a 
moderate increase in V leads mainly to a transfer of spec- 
tral weight to excitations with smaller energy loss. How- 
ever, a non-zero intersite Coulomb repulsion V is needed 
to obtain spectra related to the experiment. This fact 
led to the conclusion that the spectral intensity at the 
zone boundary is due to an exciton formation.^ On the 
other hand, the large differences between the interpreta- 
tions of the experimental spectra for Sr2Cu03 using the 
one-bandEI and the multi-band model imply that oxygen 
degrees of freedom are important for a realistic descrip- 
tion of charge excitations in the cuprates. 

In conclusion, we have carried out an investigation of 
the EELS spectrum for the one-dimensional CUO3 chain 
using an extended multi-band Hubbard model and an ex- 
tended Hubbard model. Our results for the multi-band 
model show very good agreement with experimental data 
for Sr2Cu03. In contrast to former investigations, we 
can explain the width of the main structure for all val- 
ues of momentum transfer. For the multi-band model, 
only a combination of intersite Coulomb interaction and 
charge-transfer energy is relevant for the loss function. 
Consequently, we find no evidence for enhanced inter- 
site interactions in Si^CuOa. The different explanations 
for the spectral intensity at the zone boundary using the 
one-band and the more realistic multi-band model shows 
that the oxygen degrees of freedom are important for the 
description of charge excitations. 
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